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A  Hydrogen-Helium  mixture  was  chosen  to  investigate  the  structure  of  a  counterflow  diffusion 
flame.  Reacting  and  non  reacting  conditions  were  studied  at  the  same  Reynolds  number.  To  study  the 
reaction  zone  structure,  high  speed  tomography  based  on  Mie  scattering  was  employed  using  a  copper 
vapor  laser  and  a  Fastax  highspeed  camera.  LDV  measurements  were  also  obtained.  Different  seeding 
techniques  were  used  to  visualize  both  the  turbulent  air  and  fuel  jets.  The  tomographic  records  were  digi¬ 
tized  and  recorded  in  a  digital  computer  for  statistical  treatment  Significant  differences  in  the  wrinkle 
scales  between  the  reacting  and  non  reacting  flows  were  found.  A  fractal  statistical  analysis  of  the  tomog¬ 
raphy  records  was  done  to  quantify  these  differences.  Seeding  of  both  fuel  and  air  jets  provided  a  mean  for 
the  evaluation  of  the  reaction  zone  thickness.  The  strain  of  the  reaction  zone  was  obtained  from  the  time 
resolved  tomographic  records.  Local  flame  extinction  and  reignition  were  observed  for  different 
H2/He!ium  fuel  mixtures.  ^  ,  / ,  , 
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Historically,  the  evaluation  of  mixture  fraction  in  combusting  flow  has  required  multiple  concentra¬ 
tion  measurements.  However,  we  have  found  it  possible  to  choose  a  special  fuel  (85$>  He,  15%  H2)  which 
permits  the  measurement  of  the  mixture  fraction  using  only  Rayleigh  scattering.  The  Rayleigh 
intensity/mixture  fraction  correspondence  was  verified  experimentally  in  a  laminar  coflow  flame.  The 
influences  of  strain  rate  and  differential  diffusion  effects  were  investigated  by  comparing  computed  and 
measured  profiles  of  a  laminar  counterflow  diffusion  flame.  The  results  obtained  from  these  comparisons 
were  encouraging  and  suggest  that  the  established  Rayleigh/mixture  fraction  correspondence  is  valid  under 
both  turbulent  mixing  and  laminar  strained  flamelet  combustion  regimes. 


FINAL  REPORT 


AFOSR  Grant  88-0011 
Opposed  Jet  Turbulent  Diffusion  Flames 
Principal  Investigator:  L.  Talbot 
University  of  California  at  Berkeley 


Objective:  The  objective  of  this  research  program  has  been  to  investigate  experimentally  the 
dynamics  and  structure  of  opposed  jet  turbulent  diffusion  flames,  where  a  jet  of  fuel  is 
directed  counter  to  an  oncoming  airstream.  This  configuration  is  of  fundamental  interest  in 
exploring  the  effects  of  strain  and  heat  release  on  the  fluid  mechanics  of  turbulent  diffusion 
flames,  and  additionally  exhibits  some  of  the  characteristics  of  boundary  layer  separation 
caused  by  transverse  fuel  injection  in  certain  types  of  combustors. 

Principal  results:  The  majority  of  the  experiments  were  carried  out  using  a  fuel  jet  composed 
of  a  mixture  of  hydrogen  and  helium,  whose  composition  could  be  varied,  while  maintaining 
essentially  constant  jet  Reynolds  number,  to  explore  the  effects  of  burning  rate  and  heat 
release  on  the  flame  structure. 

The  experimental  techniques  employed  included  high  speed  laser  tomography  using  a 
copper  vapor  laser  synchronized  with  a  Fastax  camera  with  silicone  oil  aerosol  and  titanium 
dioxide  particles  as  the  scattering  media,  laser  velocimetry,  and  Rayleigh  scattering.  The 
scattering  media  permitted  the  imaging  separately  of  the  fuel  stream,  the  airstream  and  the 
two  in  combination.  Fractal  analysis  was  performed  on  the  flame  interface,  to  determine  its 
length  scales,  and  individual  flame  structures  could  be  identified  and  followed  in  time  to 
determine  the  local  stretch  rate  of  the  flame  interface.  Joint  pdf’s  of  axial  and  transverse 
velocity  components  were  constructed  which  provided  information  on  entrainment  rates. 

The  fractal  analysis  revealed  that  the  smallest  scales  or  wrinkle  dimensions  of  the  flame 
front  are  increased  by  heat  release,  as  are  the  largest  scales.  In  the  stagnation  region  of  the 
flame,  the  flame  front  thickness  is  comparable  to  that  of  a  laminar  diffusion  flame  and  corre¬ 
lates  well  with  the  length  scale  associated  with  the  local  strain  rate,  which  is  largest  in  the 
stagnation  region.  In  contrast,  the  entrainment  rate  away  from  the  stagnation  region  correlates 
with  the  largest  scales  associated  with  the  vortex  structure  of  the  shear  layer.  The  strain  rate, 
which  is  highest  in  the  stagnation  zone,  increases  with  heat  release  and  a  correlation  exists 
between  strain  rate  fluctuations  and  fluctuations  in  reaction  zone  position. 

At  the  lowest  concentrations  of  hydrogen  in  the  hydrogen-helium  fuel  mixture  it  was 
possible  to  produce  local  extinction  at  the  stagnation  point  of  the  flame,  even  though  combus¬ 
tion  was  maintained  in  the  shear  layer  downstream  of  the  stagnation  zone.  This  was  even 
more  dramatically  exhibited  in  the  case  of  methane-air  flames  which  have  a  lower  extinction 
strain  rate,  and  for  which  a  large  non-reacting  "hole"  in  the  flame  was  produced  at  high 
strains.  It  is  evident  that  the  opposed  jet  turbulent  diffusion  flame  is  an  excellent 
configuration  for  the  investigation  of  strain-induced  extinction  phenomena. 

Parallel  with  our  study  of  the  opposed  jet  flame,  we  investigated  the  possibility  that  a 
particular  fuel  could  be  identified  such  that  its  Rayleigh  scattering  cross-section  was  an  unique 
monotonic  function  of  mixture  fraction.  We  found  that  a  15%  hydrogen  -  85%  helium  fuel 
mixture  meets  this  criterion  and  we  demonstrated  by  means  of  measurements  in  a  laminar 
opposed  jet  (Tsugi)  flame  that  calculated  and  measured  mixture  fractions  were  in  excellent 


agreement.  This  opens  up  wide  possibilities  for  mixture  fraction  related  measurements  in  tur¬ 
bulent  diffusion  flames  using  this  fuel. 

Personnel. 

In  addition  to  the  Principal  Investigator,  the  following  personnel  were  involved  in  the 
program. 


Dr.  Philippe  J.  Goix  -  postdoctoral  researcher 
Mr.  Keith  R.  Leonard  -  graduate  student 

Dr.  Robert  K.  Cheng  -  Lawrence  Berkeley  Laboratory  staff  scientist 
Dr.  Ian  G.  Shepherd  -  Lawrence  Berkeley  Laboratory  staff  scientist 

Publications  and  Presentations 


Two  papers,  based  on  the  work  described,  are  being  prepared  for  publication.  They  are: 

"Experimental  Study  of  Turbulent  Counterflow  Diffusion  Flame  Structure",  by  P.J.  Goix,  K.R. 
Leonard  and  L.  Talbot. 

"Direct  Measurement  and  Mixture  Fraction  in  Reacting  Flow  Using  Rayleigh  Scattering",  by 
P.J.  Goix,  K.R.  Leonard,  L.  Talbot  and  J.Y.  Chen. 

Additionally,  an  invited  presentation  of  earlier  flame-stabilization  work  was  given  at  the 
ICASE  Workshop  on  Turbulent  Combustion,  October  1989,  the  proceedings  of  which  will  be 
published  by  Springer-Verlag: 

"Flameholding  in  Unconfined  Turbulent  Premixed  Flames",  by  I.G.  Shepherd,  J.R.  Hertzberg 
and  L.  Talbot. 

Poster  presentations  of  the  work  on  counterflow  turbulent  diffusion  flames  and  mixture  frac¬ 
tion  measurements  were  given  at  the  recent  23rd  Symposium  (International)  on  Combustion  in 
Orleans,  France,  August  1990. 

Preprints  of  the  above  mentioned  papers  and  abstracts  of  the  poster  presentation  are 
attached. 

Significant  Technical  Interactions 

Collaborations  with  Drs.  J.Y.  Chen  and  R.  Dibble  of  the  Sandia  Combustion  Research 
Facility  were  established  and  are  still  maintained.  Significant  technical  interactions  with  Pro¬ 
fessor  R.  Bilger  of  the  University  of  Sydney  occurred  during  his  semester-long  visit  at  Berke¬ 
ley.  We  also  participated  by  invitation  in  the  ICASE  Combustion  Workshop. 

Inventions 


None. 
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Experimental  Study  of  Turbulent  Counterflow  Diffusion 

Flame  Structure 


P.J.  Goix,  K.R.  Leonard  and  L.  Talbot 

Mechanical  Engineering  Department, 
U.C.  Berkeley,  CA  94720. 


ABSTRACT 


Laminar  counterflow  combustion  configurations  have  been  studied  extensively  for  the  investigation  of 
fundamental  combustion  phenomena,  such  as  fuel  dilution,  strain  rate,  and  pressure  dependences  on  flame 
extinction  limits.  However  very  few  studies  have  been  performed  on  turbulent  counterflow  diffusion  flames. 
The  turbulent  opposed  flow  configuration  chosen  here  consists  of  a  fuel  jet  pipe  of  d=6.4mm  in  diameter 
directed  downward  toward  an  air  jet  50  mm  in  diameter,  shielded  with  a  concentric  air  jet  of  130mm 
diameter.  The  vertical  distance  between  the  fuel  jet  exit  and  the  air  flow  can  be  adjusted  to  investigate  the 
influence  of  the  overall  strain  rate  along  the  centerline  of  this  configuration.  A  Hydrogen-Helium  mixture 
has  been  chosen  to  investigate  the  dilution  effects  on  the  structure  of  the  diffusion  flame.  To  investigate  the 
combustion  effects  on  the  turbulent  mixing,  reacting  and  non  reacting  conditions  have  been  studied  at  the 
same  Reynolds  number.  In  order  to  study  the  reaction  tone  structure,  high  speed  tomography  based  on 
Mie  scattering  has  been  employed  using  a  copper  vapor  laser  and  a  Fastax  high  speed  camera.  Different 
seeding  techniques  have  been  used  to  visualize  both  the  turbulent  air  and  fuel  jets.  The  tomographic 
records  were  digitized  and  recorded  in  a  digital  computer  for  statistical  treatment.  Significant  differences 
in  the  wrinkle  scales  between  the  reacting  and  the  non  reacting  flows  were  found.  A  fractal  statistical 
analysis  of  the  tomography  records  has  been  done  to  quantify  these  differences.  Seeding  of  both  fuel  and 
air  jets  provided  a  mean  for  the  evaluation  of  the  reaction  zone  thickness.  From  the  time  resolved  tomo¬ 
graphic  records  an  evaluation  of  the  stretch  of  the  reaction  zone  boundaries  has  been  performed. 

INTRODUCTION 

An  opposed  jet  turbulent  diffusion  flame  is  an  interesting  combustion  configuration  for  several  rea¬ 
sons.  One  is  that  it  facilitates  the  study  of  the  effect  of  strain  on  the  reaction  zone  in  tb  stagnation  region 
and  thus  extinction  phenomena  can  be  readily  investigated  even  under  conditions  when  co-flowing  tur¬ 
bulent  diffusion  flames  may  be  unstable.  A  second  interesting  feature  of  this  configuration  is  that  it  incor¬ 
porates  at  least  qualitatively  some  of  the  characteristics  of  the  shock  wave  induced  boundary  layer  separa¬ 
tion  which  occurs  ahead  of  transversally  injected  fuel  jets  such  as  those  which  have  been  considered  for 
employment  in  scramjet  combustors  and  indeed,  the  opposed  jet  configuration  has  itself  been  mentioned  as 
a  possible  candidate  for  scramjet  fuel  injection  although  there  are  obvious  drag  penalties  involved.  Addi¬ 
tionally,  the  configuration  provides  a  challenging  problem  in  the  evolution  of  a  combustion  zone  from  a 
turbulent  stagnation  combustion  zone  to  an  entrainement  controlled  shear  layer  region.  Laminar  one¬ 
dimensional  stagnation  diffusion  flames  have  been  studied  theoretically  and  experimentally  quite  exten¬ 
sively  recently  [l-3j.  In  particular  the  dilution  effect  and  strain  effect  on  the  heat  release  variations  and 
extinction  limit  conditions  have  been  emphasized  (2-3).  However  very  few  studies  have  been  done  on  such 
flames  in  turbulent  flow  conditions.  How  the  turbulent  strain  rate  at,  and  the  bulk  strain  rate  a»,  can  affect 
the  reaction  zone  structure  of  the  flames  is  not  well  known,  and  the  effect  of  turbulence  on  the  mixing 
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process  in  a  stagnation  flame,  is  still  an  important  issue. 

Our  purpose  is  to  address  this  problem  by  studying  experimentally  a  two  opposed  jets  configuration  in 
reacting  and  non-reacting  conditions.  In  order  to  analyze  the  two  boundaries  of  the  reaction  zone  a  high 
speed  tomography  laser  diagnostic  has  been  used  to  investigate  an  hydrogen/air  turbulent  diffusion  flame. 
As  hydrogen  flame  burns  very  fast,  helium  dilution  has  been  used  to  slow  down  the  chemistry.  The  advan¬ 
tages  of  choosing  He  for  dilution  are  first  to  minimize  the  preferential  diffusion  effect  between  the  hydrogen 
and  the  diluent  and  second  to  keep  the  dynamic  flow  field  conditions  as  constant  as  possible  when  helium 
dilution  effects  are  to  be  studied.  Reynolds  number  effects  and  momentum  fluxes  remaining  basically  con¬ 
stant  for  different  H2/He  mixtures  at  constant  mass  flow  rates.  In  order  to  investigate  the  dynamics  of 
such  a  turbulent  flowfieid  two  component  LDV  measurements  were  carried  out  in  reacting  and  non- 
reacting  conditions  for  the  same  Reynolds  number. 

EXPERIMENTAL  DETAILS 

The  Opposed  Jet  Burner 

The  combustion  configuration  investigated  consisted  of  a  fuel  jet  pipe  of  d=6.4mm  diameter  directed 
downward  toward  an  air  jet  of  50mm,  shielded  by  130mm  air  jet  coflow.  In  order  to  prevent  preheating  of 
the  fuel  at  the  jet  exit,  a  small  diameter  cooling  coil  (d=2mm)  was  wrapped  around  the  fuel  pipe  and,  the 
whole  fuel  pipe  was  shielded  with  a  ceramic  tube.  The  air  jet  and  the  air  coflow  velocities  have  been 
matched  at  5  m/s.  Turbulence  (u’/U=  8%)  in  the  air  flow  was  generated  by  a  perforated  plate.  The  verti¬ 
cal  distance  between  the  fuel  jet  exit  and  the  air  flow  could  be  adjusted  to  vary  the  overall  strain  rate  in 
the  reaction  zone.  Hydrogen  fuel  with  various  levels  of  helium  concentration  were  used  to  study  the  effect 
of  dilution  on  the  extinction  of  the  counterflow  diffusion  flame.  To  investigate  the  structure  and  the 
dynamics  of  the  flame,  high  speed  tomography  and  laser  Doppler  velocimetry  were  used  and  will  be 
described  in  the  following  section. 

High  Speed  Laser  Tomography 

The  laser  tomography  technique  is  based  on  Mie  scattering  and  has  been  used  extensively  over  the  last  ten 
years  [4-8].  A  laser  light  sheet  of  60mm  width  and  .5mm  thickness  is  produced  by  a  copper  vapor  laser  and 
two  cylindrical  lenses  as  shown  as  Fig  (1).  The  copper  vapor  laser  provides  5  mJ  pulses  of  20ns  duration  at 
adjustable  repetition  rates  up  to  6khz  at  X=5Ilnm.  To  follow  the  high  repetition  rate  of  the  laser  we  used 
a  high  speed  Fastax  camera  to  image  the  laser  sheet.  The  Fastax  image  synchronization  pulses  are  used  to 
trigger  the  laser.  The  laser  is  fired  every  two  camera  pubes,  to  record  the  radiation  illumination  back¬ 
ground  coming  from  the  flame.  The  film  used  to  record  the  tomography  events  was  Kodak  4X  reversal  (400 
ASA).  Two  types  of  seeding  particle  were  used,  titanium  dioxide  (Ti02)  and  silicone  oil  particles.  Titanium 
dioxide  particles  were  produced  by  the  reaction  TiCU  +  2H20  —  >  Ti02  +  4(HC1).  By  seeding  the  fuel  jet 
with  TiCl4  and  H20  a  very  dense  cloud  of  Ti02  can  be  obtained.  This  technique  has  been  used  by  various 
authors.  Zur  Loye  and  Bracco  [7],  investigated  the  effects  of  the  presence  Ti02  and  HC1  on  the  combustion 
process  and  found  it  to  be  negligible.  Silicone  oil  and  (TiCU  +  H20)  were  injected  respectively  in  the  air 
stream  and  the  fuel  stream  using  an  atomizer  through  a  bypass  circuit.  To  determine  at  which  tempera¬ 
ture  the  maximum  Mie  scattering  gradient  corresponds  for  each  type  of  particles,  a  cylindrical  laminar 
diffusion  flame  burner  was  set  up.  The  fuel  and  the  air  jets  were  seeded  respectively  with  Ti02  and  silicon 
oil  and  a  Mie  scattering  point  experiment  similar  to  Kennedy  and  Kent  [9]  was  carried  out.  Point  tem¬ 
perature  and  Mie  scattering  measurements  were  recorded  at  the  same  locations  using  a  Pt/PtlQ%Rh  ther¬ 
mocouple  of  300  fitn  bead  mean  diameter.  Fig(2)  shows  the  comparison  of  the  two  normalized  Mie  scatter¬ 
ing  profiles  and  the  temperature  profile  taken  at  10mm  from  the  burner  exit.  One  can  note  the  very  sharp 
intensity  decrease  within  the  laminar  reaction  zone,  since  for  tomography  visualization  our  concern  is  the 
maximum  intensity  gradient.  The  temperatures  corresponding  to  this  criterion  for  both  particle  types  are 
found  to  be  480K  for  the  silicone  oil  and  420K  for  Ti02.  A  high  speed  tomography  negative  film  frame  is 
projected  and  is  imaged  by  a  video  camera.  The  video  signal  is  digitized  and  stored  on  a  512X512  I2S 
matrix  memory  with  258  grey  levels.  The  scaling  factors  in  the  X  and  Y  directions  are  .125  mm  per  pixel 
and  and  .105  mm  per  pixel  respectively.  After  background  removal  a  threshold  corresponding  to  the  max¬ 
imum  intensity  gradient  is  determined  from  the  grey  level  probability  density  function  and  the  edge 
corresponding  to  this  threshold  is  computed  and  stored  for  further  statistical  treatment. 
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Velocity  Measurements 

To  measure  the  average  flow  field  of  the  turbulent  counterflow  H2-diluted  Helium/Air  flame  two  com¬ 
ponent  LDV  point  measurements  were  performed.  Details  of  the  two  LDV  components  and  the  data 
acquisition  are  included  in  Cheng  [9].  A  Spectra-Physics  4  watt  argon-ion  laser  was  the  light  source  for  the 
two  component  (488nm  and  514nm)  LDA  system  which  uses  two  TSI  model  1980A  frequency  counters. 
Aluminium  dioxide  particles  of  .3  pm  in  diameter  were  chosen  to  track  both  the  fuel  and  the  air  jets. 
Since  the  dynamic  range  of  the  velocity  field  in  the  counterflow  diffusion  flame  was  expected  to  be  between 
60m/s  and  0  for  the  fuel  jet  region  and  in  between  0  and  12m/s  in  the  product  and  cold  air  region  a 
differential  frequency  shifting  of  5  Mhz  for  both  components  was  chosen.  To  reduce  the  velocity  signal  fre¬ 
quency  and  the  shifting  requirements  the  two  LDV  channel’s  optics  set  up  were  rotated  45  degrees  from  the 
burner  axis.  The  two  velocity  components  were  recored  simultaneously  and  transformed  to  the  original 
burner  frame  of  reference  for  further  statistics.  The  number  of  samples  for  each  records  was  8192. 

RESULTS 

The  Two  Opposed  Jet  Velocity  Field 

The  turbulent  counterflow  configuration  was  studied  under  non  reacting  and  reacting  conditions.  In  both 
cases  the  fuel  jet  and  air  jet  exits  was  separated  by  190mm.  The  hydrogen  fuel  jet  was  diluted  with  60%  of 
Helium  in  volume.  The  Reynolds  numbers  based  on  the  maximum  jet  velocities  were  2700  and  3300  respec¬ 
tively  for  the  non-reacting  (N.R.)  and  reacting  (R.)  flow  field.  Figs  3  a  ,b  represent  the  average  velocity 
vector  fields  of  the  turbulent  opposed  jets  configurations  in  non-reacting  and  reacting  conditions.  In  these 
figures  the  the  velocity  vectors  are  plotted  in  the  X,Y  plane  where  the  jet  centerline  corresponds  to  Y=0, 
and  X=0  corresponds  to  the  fuel  jet  exit.  In  both  cases  the  velocity  field  is  axisymetric.  One  can  note  the 
difference  in  the  penetration  depth  of  the  fuel  jet  between  the  N.R  (Xp=45mm)  and  the  R.  (Xp=110mm) 
cases.  The  divergence  of  the  flow  field  near  the  stagnation  region  is  much  stronger  in  combustion.  These 
differences  in  penetration  depths  and  divergence  are  mainly  associated  with  volumetric  expansion  occuring 
during  combustion.  It  is  also  clear  from  these  two  figures  that  the  air  flow  streamlines  diverge  near  the 
stagnation  region  and  then  rotate  as  they  are  entrained  in  large  scale  structures  by  the  fuel  jet.  The  size  of 
these  large  scale  structures  differ  from  the  R.  to  the  N.R  cases,  being  approximatively  40mm  and  20mm 
respectively.  In  Figs  (4  a,b)  are  plotted  the  centerline  velocity  profiles  in  both  cases.  In  both  cases  the  fuel 
and  the  air  jets  have  been  seeded  separately  with  A1203  particles.  One  can  note  significant  differences 
between  the  two  cases  even  though  both  the  air  jet  exit  and  the  fuel  jet  exit  velocities  are  matched  to 
+5m/s  and  -65m/s  respectively.  The  average  strain  rates  deduced  from  these  profiles  at  different  jet  loca¬ 
tions  are  given  in  Table  I.  At  the  stagnation  region  the  average  strain  rate  values  are  found  to  be  800«~‘ 
(fuel  side)  and  55«_1  (air  side)  for  the  N.R  case  and  500  «~l  reacting  (fuel  side)  and  95a-1  (air  side).  In 
order  to  have  a  better  understanding  of  the  flow  field,  joint  probability  density  functions  of  the  two  velo¬ 
city  components  have  been  computed  in  both  the  shear  and  stagnation  region  in  reacting  conditions.  In 
each  case  the  fuel  and  the  air  jets  have  been  seeded  separately.  Figs  5  a,b  represent  joint  pdf  contours, 
jpdf(u,v)  computed  at  the  stagnation  region.  These  contours  correspond  to  70%,  40%  and  20%  of  the 
maximum  peak  value.  For  each  case  Umean  and  Vmean  axes  are  plotted.  A  comparison  of  Fig  5a  and  5b 
show  a  dear  bimodal  behavior  around  the  zero  velocity  due  to  the  intermittency  of  the  stagnation  region. 
Since  the  velocity  gradient  is  much  steeper  in  the  fuel  side,  the  jpdf  is  much  broader.  Due  to  the  intermit¬ 
tency  even  stationary  particles  in  the  stagnation  region  can  give  positive  or  negative  velocity  values,  and 
one  has  to  be  careful  in  the  determination  of  the  average  strain  rate  at  the  stagnation  region.  Figs  6  a-e 
and  Figs  7  are  represent  the  jpdf(u,v)  contours  measured  and  computed  in  the  shear  region  where  the 
combustion  is  taking  place  at  -90mm  and  -100mm  of  the  fuel  jet  exit  (left  side  only,  negative  Y).  The  posi¬ 
tions  at  which  the  jpdf(u,v)  have  been  calculated  are  indicated  on  the  velocity  plots  with  a  circle  (Fig  3.a). 
The  negative  and  positive  transverse  velocities  point  toward  the  outside  and  inside  the  fuel  jet  respectively 
as  shown  on  Fig  8  where  the  four  quadrants  and  their  meanings  are  designated.  In  each  case  the  Umean 
and  Vmean  axes  are  plotted.  Figs  6  a-e  indicate  that  both  fuel  and  air  fluid  are  rotated  from  quadrant  1  to 
quadrant  2,  from  X— 90m  Y«-15mm  to  X—-90mm  Y— 9mm,  entraining  hot  products  in  the  fuel  side.  Figs 
7  a-e  represent  the  same  jpdf(u,v)  for  the  same  Y  locations  but  in  this  case  at  X~-100mm.  The  joint  pdf 
velocity  contours  are  narrower  particularly  for  the  air  seeded  cases,  and  the  rotation  phenomenon  pointed 
out  previously  is  less  pronounced.  However  we  can  still  observe  a  20%  peak  contour  in  the  negative 
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longitudinal  velocity  quadrant,  Fig  7  (a),  indicating  the  beginning  of  the  rotation  from  quadrantl  to  qua- 
drant2  representative  of  entrainment. 

Tomography  Results 

Figs  9  a,b  show  instantaneous  tomography  snapshots,  where  the  fuel  only  has  been  seeded  with  Ti02  par¬ 
ticles.  In  both  cases  the  images  are  bimodal  in  grey  levels  and  the  transition  between  the  black  region  (air) 
and  the  bright  one  (fuel)  is  very  sharp.  In  the  cold  condition  this  transition  corresponds  to  the  fuel/air  jet 
interface  and  in  combustion  this  interface  is  assumed  to  be  an  isotherm  interface  at  Ti=420K  as  shown  in 
Fig  2.  These  interfaces  have  been  recorded  in  a  PDPll  computer  memory  using  a  threshold  continuity 
algorithm.  In  order  to  understand  the  wrinkle  process  of  the  highly  convoluted  interfaces,  a  fractal  algo¬ 
rithm  has  been  applied  to  the  tomographic  records.  Fractal  analysis  has  been  recently  used  in  turbulent 
combustion  [11-15).  It  has  been  shown  by  Mandelbrot  jll)  that  a  highly  convoluted  fractal  surface  can  be 
described  by  three  parameters  which  are  the  fractal  dimension  D,  the  inner  cutoff  £,  and  the  outer  cutoff 
£„.  These  three  parameters  can  be  interpreted  respectively  as  the  wrinkle  scale  hierarchy  of  the  surface,  the 
smallest  wrinkle  scale  and  the  largest  wrinkle  scale  of  the  surface.  The  fractal  algorithm  used  here  is  a 
counter  segment  algorithm.  This  algorithm  is  similar  to  those  used  by  Mandelbrot  [ll]  and  Gouldin  et  al 
[12]  to  analyze  flame  fronts  in  premixed  combustion  and  is  described  in  [13].  Basically  the  length  of  the 
digitized  curve,  repiesenting  either  an  isotherm  or  the  fuel/air  jet  interface,  is  measured  with  a  straight 
line  segment  c  which  varies  in  size.  The  length  L(«)  is  obtained  by  the  relationship  L(c)=N(e)(  To  minimize 
noise  problems  statistics  over  30  samples  have  been  processed.  Figs  10  a,b  show  a  comparison  of  the  fractal 
plots  corresponding  to  the  typical  tomography  contours  shown  Figs  9  a,b.  Both  of  the  plots  exhibit  a  frac¬ 
tal  character,  the  evolution  of  the  length  versus  the  line  segment  size  in  log-log  space  is  linear  in  both  cases 
and  the  fractal  dimension  is  found  to  be  equal  to  D*=2.18  (+/-.005)  for  the  N.R.  case  and  D=2.07  (+/-.005) 
for  the  R.  case.  The  plots  are  linear  for  at  least  one  decade  of  the  log  space.  The  inner  cutoff  and  the  outer 
cutoff  are  determined  where  the  curves  depart  respectively  at  small  and  large  t  from  the  fractal  slope.  For 
the  N.R.  case  £,•  =  1.3mm  and  £„  =  16mm  while  for  the  R.  case  «,  =  2mm  and  („  —  20  mm.  One  can  note 
significant  differences  between  the  N.R  and  the  R.  stagnation  regions.  Overall  the  fractal  dimension  is 
found  to  be  much  smaller  in  the  reacting  case,  indicating  a  smoothing  effect  on  the  isotherms  of  the  flow 
field  due  to  the  combustion.  Another  interesting  feature  to  note  is  that  both  the  inner  cutoff  and  the  outer 
cutoff  in  the  reacting  flow  are  larger  than  the  non  reacting  case  which  indicates  that  the  turbulent  velocity 
spectra  in  combustion  are  shifted  toward  larger  scales. 

Figs  11  a,b  represent  two  instantaneous  tomography  snapshots  In  this  case  both  the  fuel  and  air  jets 
have  been  seeded  with  Ti02  and  silicone  oil  particles.  The  reaction  zone  (black  region)  is  clearly  delineated 
by  the  two  isotherms  T,„2,  T^.  The  determination  of  the  position  of  these  two  isotherm  contours  leads  to 
a  spatial  determination  of  the  reaction  zone  thickness.  Qualitatively,  the  tomography  pictures  show  a 
significant  increase  of  the  reaction  thickness  from  the  stagnation  region  toward  the  shear  region.  One  can 
also  see  the  vortex  development  along  the  shear  region  and  entrainment  of  hot  product  gas  in  the  fuel  jet. 
From  time  sequences  typical  of  Fig  11  a,b  it  is  possible  to  isolate  a  turbulent  structure  in  both  the  shear 
and  stagnation  regions  and  to  determine  their  length  evolutions  with  time.  However  the  criterion  chosen  to 
isolate  a  turbulent  structure  is  subjective,  and  we  decided  to  isolate  a  turbulent  structure  when  it  is  bound 
by  a  singularity  (second  derivative  discontinuity)  at  each  end.  Following  the  singularities,  we  are  then  able 
to  measure  their  lengths.  This  leads  to  a  one  dimensional  evalution  of  the  stretch  factor 
5<r,«l/L(t)  *  dL(t)/dt,  and  assuming  isotropic  stretching  of  the  surface  area  the  stretch  factor  can  be 
written  as  5<r2=l/A(t)  *  dA(t)/dt=2  *  Strx .  The  results  obtained  from  the  high  speed  tomography  pic¬ 
tures  for  the  stagnation  region  and  the  shear  zone  will  be  presented  separately. 

Stagnation  Zone 

The  determination  of  stretch  and  the  instantaneous  laminar  reaction  zone  thickness  evaluation  is  very 
interactive  and  quite  laborious,  and  for  this  reason  at  present  statistics  of  those  quantities  over  only  30 
independent  samples  are  presented.  One  sample  correspond  to  two  snapshots  .57ms  apart.  The  average 
turbulent  reaction  zone  thickness  is  found  to  be  e„«1.5mm  with  an  rms  of  .5  mm.  The  average  stretch 
rate  S(r2  deduced  is  equal  to  1500s"-1  with  an  rms  of  650  »~l.  This  rms  could  probabily  be  reduced  by 
increasing  the  number  of  samples  in  the  statistics.  For  comparison,  the  stretch  rate  has  been  also  computed 
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at  the  stagnation  region,  for  the  same  number  of  samples, in  the  non  reacting  flow  case  (N.R).  The  average 
stretch  rate  value  obtained  in  the  non  reacting  case  is  Str2  =  820  s~l  with  an  rms  of  450  a-1.  This  value  is 
significantly  lower  than  the  stretch  rate  in  the  N.R  conditions.  These  differences  will  be  discussed  later.  It 
is  however  worth  noting  that  the  strain  rate  determination  requires  a  highly  convoluted  contour 
(representing  the  highest  strain  rate)  and  the  statistical  evaluation  of  the  stretch  may  correspond  to  the 
upper  bound  of  the  turbulent  stretch  distribution.  Fig  12  a  represents  the  pdf  of  the  position  of  the  lean 
side  reaction  zone  (air  boundary),  marked  by  the  isotherm  T*  =480 K.  This  pdf  has  been  obtained  from 
200  time-independent  samples  along  the  centeiline.  From  this  pdf  the  rms  of  the  reaction  zone  position  in 
the  stagnation  zone  is  determined  to  be  3.5  mm. 

Recirculation  Zone 

Fig  12  b  represents  the  evolution  of  the  pdf  of  position  of  along  the  shear  layer.  The  rms  of  position 
deduced  from  the  pdf  decreases  from  3.25  mm  at  X=  -100  mm  to  2.5  mm  at  X=  -95  mm  and  is  constant 
up  to  X=  -75  mm.  These  results  indicate  the  stabilizing  role  played  by  the  recirculation  region  in  opposing 
the  jet’s  lateral  movements.  As  has  been  done  for  the  stagnation  region,  the  stretch  factor  has  also  been 
determined  along  the  shear  layer.  It  is  clear  Fig  lib  that  while  the  shear  turbulent  structures  are  being 
rolled  up  they  are  convected  downstream.  The  stretch  factor  during  the  roll  up  varies  between  100  «_1  and 
800  «-1  while  the  convection  speed  of  the  structure  varies  between  6  m/s  and  13  m/s  where  the  velocity 
vector  always  diverges  from  the  axis  of  symmetry. 

DISCUSSION 


The  mixture  fraction  variable,  Z,  a  conserved  scalar,  is  very  useful  in  the  modeling  of  the  internal 
structure  of  diffusion  flames  [16-18] .  By  definition,  the  mixture  fraction  is  bounded  by  0  in  pure  oxidant 
and  by  1  in  pure  fuel.  Assuming  fast  chemistry  and  thermodynamic  equilibrum,  the  temperature  of  the 
flame  in  mixture  fraction  space  can  be  obtained.  Fig  12  shows  the  temperature  profile  in  mixture  fraction 
space  for  the  conditions  of  the  H2-helium/air  diffusion  flame  studied.  The  temperature  and  mixture  frac¬ 
tion  at  the  stoichiometric  surface  are  T=2000  K  and  Z,  =  0.1  respectively.  The  problem  is  to  find  the 
correspondance  between  the  mixture  fraction  and  physical  space.  As  shown  by  Bilger  [16],  Peters  [17]  and 
Williams  [18]  ,  using  the  fuel  and  oxidant  mass  fraction  conservation  equations  one  can  derive  an  equation 
for  the  mixture  fraction  under  a  constant  diffusivity  approximation, 

P^.  +  pvVZ=V(pDVZ)  (1) 


For  a  laboratory  coordinate  system  in  the  case  of  a  steady  laminar  stagnation  flame  this  equation  can  be 
written  as  follows: 


dZ  d  ,  ndZ, 

fuT,  m  -  T.^T^ 


(2) 


Assuming  pD  constant  and  a  constant  strain  rate  a,  one  can  obtain  an  analytical  solution  for  the  mixture 
fraction: 


*(*)- 


(3) 


where:  a 


a.  1 

2D’  V7 


has  a  length  scale  dimension  which  can  be  related  to  the  reaction  zone  thickness  of 


the  laminar  diffusion  flame.  Our  purpose  here  is  to  compare  the  instantaneous  thickness  of  the  reaction 
zone  of  the  turbulent  opposed  jet  diffusion  flame  determined  by  high  speed  tomography,  with  the  thickness 
calculated  for  the  laminar  case  under  the  same  strain  rate  conditions.  The  comparison  shown  in  Fig  2  of 
the  maximum  intensity  gradient  with  a  temperature  profile  in  a  laminar  case  for  Tio2  and  oil  particles 
provides  a  means  of  marking  the  isomixture  fraction  values  and  Zrf.  Using  the  temperature/mixture 
fraction  plot  of  Fig  13  as  a  look  up  table,  we  can  evaluate  the  isomixture  fraction  contours  corresponding 
to  the  oil  and  TI02  maximum  intensity  gradient  contours,  which  are  Z ^  «  .001  and  Zj y#J  =  .8.  Several 
experimental  and  theoretical  studies  of  laminar  opposed  jet  reacting  flows  have  shown  that  the  laminar 
velocity  profile  in  the  reaction  zone  is  not  monotonic.  In  particular,  Tsuji  [1]  and  more  recently  Pellet  et  al 
[3]  have  shown  that  the  typical  velocity  profile  increases  linearly,  crossing  the  stagnation  region,  from  the 
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jet  side  to  the  flame  (the  stochiometric  mixture  fraction)  and  then  decreases  sharply  in  a  region  called  the 
"upstream  boundary  of  flame  front"  and  finally  increases  linearly  on  the  air  jet  side.  In  all  cases  the  velo¬ 
city  gradients  on  the  air  side  are  much  smaller  than  the  gradients  on  the  fuel  side.  The  sharp  gradient 
between  the  flame  and  its  upstream  boundary  position  is  attributed  to  the  heat  release  and  usually  occurs 
in  a  very  thin  region  typically  0.2  mm  in  width.  In  our  case  this  sharp  gradient  must  occur  between  Z,=.  1 
and  Z=0.  It  was  not  possible  to  resolve  this  region  with  our  velocity  measurements  due  to  the  intermit¬ 
tent  movement  of  the  reaction  zone.  However,  as  shown  in  Table  1  the  gradients  on  the  air  side  are  about 
five  time  smaller  than  on  the  fuel  side.  Figs  14  a,b  show  the  laminar  reaction  zone  profiles  deduced  from 
equation  (3).  The  diffusivity  coefficient  chosen  is  D  (He-N2)  at  T=298  K  and  also  at  T=1150K  which 
corresponds  to  (Tmax+  Tmin)/2.  The  strain  rate  has  been  varied  from  500  a-1  to  1500  s~‘  which  are 
basically  the  average  strain  rate  deduced  from  the  average  fuel  centerline  velocity  profile  and  the  stretch 
rates  deduced  from  the  tomography  length  measurements.  The  mixture  fraction  gradients  deduced  from 
these  curves  from  Z(l<)2  =  .8  to  Z  =  .1  are  given  in  Table  2  together  with  the  average  mixture  fraction  gra¬ 
dient  obtained  from  the  instantaneous  tomography  snapshots.  From  these  comparisons,  the  curves 
obtained  with  D(He-N2)  evaluated  at  Tmean=1150K  seem  to  give  results  similar  to  those  measured  while 
D(He-N2)  for  298K  leads  to  a  much  steeper  mixture  fraction  gradient.  The  best  fit  seems  to  be  at  a=1000 
s~'  which  is  between  the  average  strain  rate  500  «-1  and  the  average  stretch  rate  1500  a-1.  It  is  worth 
noting  that  the  rms  of  the  reaction  zone  thickness  is  0.5  mm  and  using  the  reaction  thickness  deduced  from 
eq  (3),  this  rms  corresponds  to  a  strain  rate  of  500  a-1.  Since  the  local  stretch  rate  was  deduced  empirically 
(when  the  stagnation  zone  was  sufficiently  wrinkled)  it  might  correspond  to  the  smallest  scales  of  the  tur¬ 
bulence  spectrum  and  thus  to  the  maximum  instantaneous  strain  rates.  This  suggests  that  the  reaction 
zone  fluctuations  can  be  correlated  to  the  strain  rate  fluctuations.  A  more  consistent  statistical  analysis  of 
both  the  reaction  zone  thickness  and  the  instantaneous  stretch  rate  is  necessary  to  confirm  this  trend.  The 
effect  of  nitrogen  dilution  on  H2/air  diffusion  flame  extinction  was  investigated  by  Dixon-Lewis  and  Mis- 
saghi  [2].  They  found  in  particular  a  significant  decrease  in  peak  temperature  beginning  for  strain  rates 
between  100  a-1  and  1000  a~\  for  an  equimolar  hydrogen/ nitrogen  mixture.  In  the  case  "f  the 
hydrogen/helium  mixture  investigated  here  no  extinction  has  been  observed.  However,  the  pdf  of  reaction 
zone  position  presented  in  Fig  12,  clearly  shows  two  modes.  Since  we  have  seen  that  the  average  flow  field 
in  combustion  exhibits  a  penetration  depth  which  is  twice  the  penetration  depth  under  cold  conditions, 
heat  release  might  induce  strong  penetration  depth  fluctuation  and  might  explain  the  two  peaks  in  the 
reaction  zone  position  pdf. 

CONCLUSION 

A  turbulent  opposed  jet  H2-He  diffusion  flame  has  been  investigated  experimentally  by  high  speed 
tomography  and  two  component  laser  velocimetry.  Different  seeding  techniques  have  been  used  to  visualize 
both  sides  of  the  reaction  zone  (fuel  and  air  sides).  A  fuel  mixture  (60%  He,  40%  H2)  at  a  Reynolds 
number  close  to  3000  has  been  investigated  in  both  reacting  and  non-reacting  conditions.  From  the  velo¬ 
city  field  measurements  and  the  tomography  visualizations  two  main  regionthave  been  identified: 

-  Shear  region:  where  the  hot  product  entrainment  into  the  fuel  jet  is  significant  and  occurs  at  large  scales. 

-  Stagnation  region:  where  the  turbulence  wrinkles  the  reaction  zone  and  moves  it  back  and  forth  around 
an  average  stagnation  point  location. 

A  fractal  segment  counting  algorithm  was  applied  to  the  tomography  records  and  the  smallest,  largest  and 
fractal  dimensions  of  the  wrinkle  scales  on  the  fuel  side  of  the  stagnation  region  have  been  deduced.  A 
local  stretch  factor  has  also  been  estimated  for  non-reacting  and  reacting  conditions.  By  simultaneous 
seeding  of  both  fuel  and  air  jets,  the  instantaneous  reaction  zone  thickness  was  obtained  and  an  average 
thickness  and  rms  were  deduced  from  a  small  number  samples.  From  these  results  we  have  found: 

1)  The  size  of  the  wrinkles  in  the  fuel  jet  stagnation  zone  are  much  larger  in  the  reacting  case  than  in  the 
nonreacting  case. 

2)  The  instantaneous  stretch  rate  is  higher  at  the  stagnation  zone  than  in  the  shear  zone. 

3)  Although  the  wrinkle  scales  are  smaller  in  the  non-reacting  case,  the  stretch  value  is  higher  in  the  react¬ 
ing  case.  This  indicates  a  possible  temperature  induced  viscosity  increase  effect  on  the  wrinkling  process. 

4)  The  reaction  zone  thickness  evaluated  at  the  stagnation  region  scales  with: 


0 


where  the  coefficients  are  chosen  as: 

D,  the  He-N2  binary  diffusion  coefficient,  evaluated  at  1150K  and 
a,  the  strain  rate,  at  1000  s'1. 

This  strain  rate  is  midway  between  the  average  strain  rate  evaluated  from  the  LDV  measurements  (500 
s'1)  and  the  upper  bound  of  the  stretch  evaluated  from  the  tomographic  time  series  (1500  s'1).  The  rrr.s  of 
the  reaction  zone  thickness  seems  to  correspond  to  a  strain  rate  rms  of  500  s'1.  This  result  indicates  a 
correlation  between  fluctuations  of  the  reaction  zone  thickness  and  strain  rate  fluctuation.  Further  statis¬ 
tics  are  necessary  to  confirm  this  trend  and  are  in  progress. 
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High  speed  laser  tomography  setup. 

Comparison  of  the  Mie  scattering  intensity  profiles  (Ti02  and  silicone  oil  particles),  with  the 
temperature  profile,  obtained  at  10  mm  from  the  burner  exit  of  a  laminar  diffusion  flame. 

Average  velocity  vector  plots:  (a)  reacting  conditions,  (b)  non-reacting  conditions. 

Centerline  velocity  profiles  in  both  non-reacting  conditions:  (a)  mean  velocity  U,  (b)  rms  u’and 
v\ 

Jpdf  (u.v)  contours  obtained  in  reacting  condition  near  the  stagnation  region  of  the  two  jets  at 
107  mm:  (a)  fuel  jet  is  seeded,  (b)  air  jet  is  seeded. 

Jpdf  (u,v)  contours  obtained  in  reacting  condition  in  the  shear  region  of  the  two  jets  at  X=  -90 
mm  for  Y=  -15  mm  to  -9  mm:  (a),  (c),  (e)  fuel  jet  seeded  and  (b),  (d),  (f)  air  jet  seeded. 

Jpdf  (u,v)  contours  obtained  in  reacting  condition  in  the  shear  region  of  the  two  jets  at  X= 
-100  mm  for  Y=  -15  mm  to  -9  mm:  (a),  (c),  (e)  air  jet  seeded  and  (b),  (d),  (f)  air  seeded. 

Schematic  representing  the  four  quadrants  of  the  jpdf. 

Instantaneous  tomography  snapshots  of  the  two  opposed  jets.  In  both  cases  only  the  fuel  is 
seeded  with  Ti02  particles:  (a)  non-reacting  conditions  (b)  reacting  conditions. 

Typical  fractal  length  plot  in  log-log  space  normalized  by  the  outer  cutoff  length,  (average  of 
30  samples):  (a)  non-reacting  conditions,  (b)  reacting  conditions. 

Instantaneous  tomography  snapshots  570  pa  apart.  Fuel  seeded  with  Ti02  and  air  seeded 
with  silicone  oil:  (a)  t=0,  (b)  t=570  pa. 

Pdf  of  position  of  T— 480K  isotherm  (air  jet  seeded  with  silicone  oil)  Reacting  conditions. 
Statistics  over  200  samples:  (a)  shear  region,  (b)  stagnation  region. 

Reaction  zone  temperature  profile  in  mixture  fraction  space. 

Mixture  fraction  profiles  near  the  stagnation  zone  using  constant  strain  rate  and  constant 
diffusivity  approximations,  (a)  D(He-N2)  evaluated  at  298K,  (b)  D(He-N2)  evaluated  at  850K. 
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Abstract 


The  mixture  fraction  variable,  is  very  useful  in  describing  reaction  zone  structure  in  non-premixed 
flames.  Extinction  limits  and  turbulent  mixing  are  often  described  as  a  function  of  this  variable.  Experi¬ 
mental  evaluation  of  £  is  valuable  in  understanding  the  influence  of  turbulent  mixing  on  the  chemistry  pro¬ 
cess.  Historically,  the  evaluation  of  mixture  fraction  in  combusting  flow  has  required  multiple  concentra¬ 
tion  measurements.  Here,  a  fuel  has  been  designed  to  permit  the  measurement  of  the  mixture  fraction 
using  only  Rayleigh  scattering.  A  Rayleigh  intensity/  mixture  fraction  correspondence  has  been  obtained 
experimentally  in  a  laminar  coflow  flame.  The  influence  of  strain  rate  and  differential  diffusion  effects  has 
been  investigated  using  a  laminar  counterflow  diffusion  flame  and  shifting  equilibrium  models.  The  results 
obtained  from  these  comparisons  are  very  encouraging  and  suggest  that  the  Rayleigh/  mixture  fraction 
correspondence  established  is  valid  under  both  turbulent  mixing  and  laminar  strained  flamelet  combustion 
regimes. 

Introduction 

The  conserved  scalar  approach  has  been  established  to  be  most  useful  for  describing  the  reaction  zone 
structure  of  turbulent  non-premixed  flames.  Bilger  [1-2],  under  assumptions  of  unity  Lewis  number  and 
constant  mass  diffusivity  coefficients,  has  shown  that  turbulent  mixing  in  the  reaction  zone  can  be  modeled 
by  the  mixture  fraction.  Furthermore  it  has  been  demonstrated  by  several  authors  [2-4],  that  the  reaction 
rate,  which  is  the  relevant  quantity  to  study  in  turbulent-combustion  interaction  process,  is  a  function  of 
X,  the  mixture  fraction  dissipation  rate.  Peters  [3],  has  suggested  that  extinction  limit  conditions  of  tur¬ 
bulent  non  premixed  flames  can  be  expressed  in  terms  of  x.  In  order  to  validate  these  theories  considerable 
experimental  mixture  fraction  data  are  still  needed.  Dibble  et  al.  (6)  have  developed  a  method  to  measure 
the  mixture  fraction  in  a  coflow  turbulent  diffusion  flame.  This  point  measurement  method  involves  rather 
complicated  and  costly  laser  diagnostic  techniques  combining  Rayleigh  scattering  and  multichannel  spon¬ 
taneous  Raman  scattering.  In  this  paper  we  describe  a  very  simple  method  to  measure  mixture  fraction 
using  only  Rayleigh  scattering.  This  method  is  based  on  the  choice  of  a  fuel  providing  a  monotonic  rela¬ 
tionship  between  the  Rayleigh  scattering  signal  and  the  mixture  fraction  variable.  The  influences  on  this 
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relationship  of  differential  diffusion  effects  and  strain  effects,  in  particular  close  to  the  extinction  limit,  are 
also  examined. 

Fuel  Design 

Rayleigh  scattering  has  been  used  extensively  over  the  past  fifteen  years  as  a  diagnostic  in  turbulent 
combustion.  Under  isobaric  conditions  the  Rayleigh  signal  can  be  expressed  as  a  function  of  the  tempera¬ 
ture,  T,  the  species  mole  fractions,  Xit  and  their  scattering  cross  sections,  <r,,  where  K  is  a  constant  involv¬ 
ing  all  the  optical  and  photodetection  parameters. 

U) 

1  1-1 

Dibble  et  al  6],  have  shown  that  depending  on  the  fuel  design  the  Rayleigh  signal  across  a  combustion 

i— * 

rone  can  provide  different  types  of  information.  If,  for  example  the  average  Rayleigh  cross  section  £V,.Y, 

i-i 

remains  constant  across  the  flame  zone  then  the  Rayleigh  signal  is  inversely  proportional  to  the  tempera¬ 
ture.  This  turns  out  to  be  the  case  for  a  mixture  of  38%  methane  and  62%  hydrogen  burning  in  the  air 
On  the  other  hand,  if  the  average  Rayleigh  cross-section  divided  by  the  average  molecular  weight  is  con¬ 
stant  across  the  flame  then  the  Rayleigh  signal  is  proportional  to  the  density,  which  is  the  case  for  a  mix¬ 
ture  of  78%  hydrogen  an  22%  argon  burning  in  the  air.  Neither  of  these  mixtures  in  combustion  (due  to 
temperature  effects)  lead  to  a  mixture  fraction  measurement  using  Rayleigh  technique. 

Using  a  shifting  equilibrium  model  (S.E.M)  one  can  determine  the  temperature  and  mole  fraction 
profiles  in  the  mixture  fraction  space  for  a  given  fuel/oxidant  mixture  and  then  from  known  Rayleigh 
scattering  cross  sections  [7],  calculate  the  Rayleigh  intensity  as  a  function  of  fusing  Eq(l).  Under  combus¬ 
tion  conditions  for  most  fuels  the  Rayleigh  signal  is  generally  bimodal  in  mixture  fraction  around  the 
stoichiometric  value  where  the  temperature  is  maximum.  This  is  the  case  for  the  propane  mixture  used  by 
E.  Effelsberg  and  N.  Peters  (9).  They  used  time  series  point  measurements  and  the  Taylor  hypothesis  to 
determine  whether  the  Rayleigh  signal  corresponded  to  the  lean  or  the  rich  side  of  the  reaction  zone.  In 
this  case  the  Rayleigh  cross  section  of  propane  is  14  times  the  air  Rayleigh  cross  section  and  the  dynamic 
range  one  has  to  resolve  between  the  product  and  the  fuel  in  order  to  have  a  good  signal/noise  ratio  over 
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the  whole  mixture  fraction  range,  is  about  100.  Another  problem  concerning  propane  fuel  is  that  under 
highly  turbulent  mixing  conditions  the  flame  becomes  sooty  and  therefore  the  Rayleigh  is  severely  contam- 
inated.  The  fuel  mixture  that  we  propose  here  (85°ic  helium,  15?£  hydrogen)  has  a  smaller  mean  Rayleigh 
cross  section  than  the  combustion  products.  Then  as  shown  in  Fig.  1  the  Rayleigh  intensity  signal  across 
the  flame  in  mixture  fraction  space  is  a  monotonic  function  and  the  mixture  fraction  can  be  determined 
without  ambiguity  from  a  Rayleigh  measurement.  The  Rayleigh  signal  ratio  between  the  pure  air  side 
£  =  0  and  the  pure  fuel  side  £  =  1  is  found  to  be  22.  One  can  note  significant  differences  in  the  curve  sensi¬ 
tivity  for  the  range  of  mixture  fraction  £  *■  0  to  £,  *  ,25  (air  side)  and  £,  *  .25  to  £  —  1..  The  problem 
arising  here  is  to  be  able  to  measure  accurately  the  Rayleigh  intensity  with  a  good  signal  noise  ratio 
within  the  full  mixture  fraction  range. 

Experimental  set  up  and  Rayleigh  measurements 

In  order  to  investigate  the  Rayleigh  scattering  properties  of  the  chosen  helium-hydrogen  mixture,  a 
laminar  coflow  diffusion  flame  was  set  up.  The  burner  consisted  of  a  fuel  pipe  of  6  mm  in  diameter  sur¬ 
rounded  by  a  coflowing  air  flow  of  30mm  diameter.  The  fuel  jet  Reynolds  number  was  Re=500  The 
helium  volumetric  concentration  (S5%)  in  the  fuel  jet  was  controlled  by  a  calibrated  flow  meter.  The  Ray¬ 
leigh  measurement  set  up  is  shown  Fig  (2).  A  focused  laser  beam  from  a  copper  vapor  laser  is  directed 
across  the  flame  using  a  classical  two  lens  system.  The  beam  waist  diameter  is  3mm  over  a  length  of  10 
mm.  The  cooper  vapor  laser  is  capable  of  a  repetition  rate  of  lOkhi  with  pulse  duration  of  20ns  and  5mJ 
energy.  The  Rayleigh  light  scattered  at  90  degree  is  spatially  filtered  and  imaged  on  the  pinhole  of  a  pho¬ 
tomultiplier  using  a  50mm  f  1/1.8  lens.  The  amplified  output  of  the  photomultiplier  is  fed  through  a  box 
car  integrator.  The  box-car  is  triggered  with  appropriate  time  delay  and  set  to  integrate  over  only  the  20 
ns  duration  of  the  laser  pulses,  thus  eliminating  noise  problems  such  as  flame  and  background  radiation 
between  pulses.  The  photomultiplier  linearity  was  checked  for  the  range  of  the  operating  conditions.  For 
each  measurement  the  laser  pulse  energy  value  is  recorded  and  ratioed  to  the  Rayleigh  scattering  signal 
value.  Fig.  3  shows  a  Rayleigh  signal  radial  intensity  profile  obtained  at  10mm  from  the  coflow  burner 
exit.  Two  profiles  1mm  above  and  below  the  laser  beam  were  also  recorded,  in  order  to  evaluate  the  aver- 


age  background  intensity  level.  The  profile  of  Fig.  3  shows  clearly  the  monotonic  behavior  of  the  normal¬ 
ized  signal  I  n  *  ~  ^  from  the  air  side  to  the  fuel  side.  One  can  also  note  the  small  background 

\  *  Rmu 

contribution  on  the  signal  level.  Fig.  4  shows  the  evolution  of  the  variance  of  the  Rayleigh  intensity  signal 
versus  its  average.  The  Rayleigh  signal  probability  distribution  at  each  point  was  found  to  be  gaussian  and 
the  linear  relationship  between  Ir  and  its  variance  indicates  that  the  photon  count  rate  obeys  Poisson 
statistics.  This  result  is  important  for  future  measurement  in  turbulent  conditions,  where  the  noise  can  be 
modeled  as  a  function  of  the  average  Rayleigh  signal  intensity  value. 

Rayleigh  measurements  in  Mixture  F raction  Space 

Laminar,  one  dimensional  diffusion  flame  structure  has  been  studied  quite  extensively  recently 

(11,14).  In  particular,  the  influence  of  the  flow  field  strain  on  the  reaction  zone  structure  is  quite  well 

known  in  the  case  of  hydrogen  chemistry  (12).  The  models  employed  in  these  studies  provide  both  spatial 

and  mixture  fraction  concentration  profiles  under  different  strain  rate  conditions.  A  common  configuration 

used  to  compute  the  one  dimensional  flame  is  the  Tsuji  burner.  One  very  important  parameter  in  this 

configuration  is  the  velocity  gradient  across  the  reaction  lone  which  is  the  convection  term  of  the  problem 

When  this  convection  term  and  hence  the  strain  rate  is  increased,  a  limit  exists  beyond  which  the  flame 

cannot  be  sustained  and  it  extinguishes.  When  the  strain  rate  is  very  low  the  problem  is  a  one  dimensional 

molecular  diffusion  driven  problem.  In  this  case  (low  strain  rate),  it  is  reasonable  to  assume  that  the  flame 

structure  in  mixture  fraction  space  is  the  same  as  that  of  a  coflow  laminar  diffusion  flame,  which  is  also  a 

molecular  diffusion  driven  flame.  Fig.  5  represents  the  laminar  flame  structure  in  mixture  fraction  space  for 

the  Tsuji  burner  configuration  under  low  strain  rate  condition  a  •  10s"1  for  the  chosen  859cHe/15coH2 

fuel.  On  this  figure  are  plotted  the  mass  fractions  of  the  species  having  significant  contribution  to  IK.  From 

those  mass  fraction  profiles  and  the  Rayleigh  cross  section  values  found  in  the  literature  (5),  one  can  com- 

-tk» 

pute  the  IR’  profile  of  the  flame  as  shown  in  Fig.  6.  One  can  note  that  both^ modeled  l\  profiles  (Fig.  6 
Tsuji  burner)  and  the  one  measured  from  the  coflow  burner  Fig.  3  exhibit  the  same  features.  They  are 
both  monotonicaliy  decreasing  with  £  and  present  a  maximum  Rayleigh  intensity  gradient  on  the  air  side. 
The  main  difference  in  the  two  profiles  is  the  thickness  of  the  reaction  tone.  In  the  stagnation  case  the 
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convection  and  the  diffusion  fluxes  are  parallel  to  each  other  while  for  the  eoflow  burner  these  fluxes  are 
perpendicular  which  causes  the  reaction  zone  in  the  coflow  case  to  be  thicker  and  also  a  function  of  down¬ 
stream  distance  from  the  burner  exit.  Hence,  to  compare  the  measured  and  computed  l\  profiles  in  terms 
of  the  space 'mixture  fraction  correspondence  we  have  rescaled  the  profiles  of  Fig.  6  and  Fig  3  using  the 
maximum  /**  gradient  thickness  as  characteristic  of  the  two  profiles.  Fig.  7  represents  the  comparison  of 
the  two  profiles  thus  rescaled  using  the  maximum  gradient  thickness  for  normalization  One  can  note  very 
good  agreement  between  the  model  and  the  data.  Then,  using  the  space/mixture  fraction  correspondence 
provided  by  the  model  one  can  now  compare  the  Rayleigh  intensities  measured  and  predicted  in  mixture 
fraction  space.  Fig.  8  shows  the  excellent  agreement  between  prediction  and  experiment  under  the  low 
strain  condition  a  »  10«-1  .  It  is  also  worth  noting  that  the  region  of  £  where  the  curves  were  forced  to 
match  (maximum  l\  gradient)  corresponds  only  to  a  range  in  mixture  fraction  less  than  1  located  on  the 
air  side  of  the  reaction  zone  where  differential  diffusion  effects  are  expected  to  be  of  secondary  importance, 
as  discussed  below. 

Strain  rate  and  differential  diffusion  effects 

Figure  8  shows  clearly  the  very  good  agreement  between  the  predicted  and  measured  /' R  in  mixture 
fraction  space.  The  question  that  we  next  address  is  to  what  extent  strain  rate  may  affect  this  agreement 
Fig.  9  shows  the  reaction  zone  structure  of  the  Tsuji  type  flame  for  our  hydrogen/helium  fuel  mixture 
computed  very  close  to  the  extinction  limit  at  a  strain  rate  of  a  ■  190«~*.  The  maximum  temperature  is 
decreased  by  about  300K  in  comparison  with  the  low  strain  rate  case.  From  these  profiles  (  a  =  1 90s — 1  ), 
one  can  compute  the  corresponding  /* *  signal  and  compare  it  to  the  low  strain  rate  condition  as  shown 
Fig.  10  Once  again  even  if  the  reaction  zone  concentration  profiles  are  quite  sensitive  to  the  strain  rate 
value,  /* R  in  the  mixture  fraction  space  remains  the  same.  This  result  is  very  encouraging  and  indicates 
that  the  Rayleigh/mixture  fraction  correspondence  determined  earlier  is  practically  independent  of  the 
strain  rate,  and  thus  it  should  be  possible  to  measure  the  mixture  fraction  across  the  reaction  zone  of  a 
laminar  strained  flamelet.  Bilger  has  shown  that  mixture  fraction  variable  is  very  well  suited  to  the  study 
of  the  turbulent  mixing  occurring  in  the  reaction  zone.  When  the  mining  is  highly  turbulent,  molecular 
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effects  are  found  to  be  negligible  and  the  hypotheses  of  unity  Lewis  number  and  constant  diffusivity 
coefficients  are  more  acceptable.  Laboratory  experiments  are  very  often  made  at  moderate  Reynolds 
number  and  molecular  transport  effects  can  have  significant  contributions.  For  our  fuel  the  mass 
diffusivities  of  nitrogen,  oxygen  and  hydrogen,  helium  are  very  different  and  differential  diffusion  effects  are 
expected  to  be  significant  15:.  Figure  11  shows  for  comparison  the  /*#  variation  with  mixture  fraction 
using  the  Shifting  Equilibrium  Model  (S.E.M)  and  the  laminar  one  dimensional  model  (where  differential 
diffusion  effects  are  included).  Here  again  the  differences  between  the  two  models  are  quite  small.  Mixture 
fraction  measurements  using  Rayleigh  scattering  for  our  chosen  Hydrogen/Helium  mixture  are  very 
promising,  particularly  for  turbulent  flows  where  the  effects  associated  with  differentia]  diffusion  will  be  of 
lesser  importance  compared  with  turbulent  transport  which  can  only  be  approximated  in  currently  avail¬ 
able  turbulence  models.  The  guidance  provided  by  measurements  made  with  our  particular  fuel  could  be 
used  to  improve  turbulence  modeling  of  combustion  involving  more  practical  fuels. 

Summary 

A  fuel  consisting  of  85%  helium  and  15%  hydrogen  has  been  designed  to  make  possible  the  experi¬ 
mental  determination  of  the  mixture  fraction.  For  this  specific  fuel  the  Rayleigh  intensity  signal  is  a  mono¬ 
tonic  function  of  the  mixture  fraction.  Rayleigh  measurements  have  been  carried  out  on  a  laminar  coflow 
burner  using  this  fuel.  The  results  show  that  the  Rayleigh  intensity  can  be  measured  with  a  very  good  sig¬ 
nal  noise  ratio  for  the  entire  mixture  fraction  range.  A  laminar  counterflow  flame  model  has  been  used  to 
transform  the  spatial  measurements  in  mixture  fraction  space  and  also  to  study  the  influence  of  strain  rate 
on  the  Rayleigh/mixture  fraction  correspondence,  (particularly,  close  to  extinction),  which  has  been  found 
to  be  negligible.  A  comparison  of  the  laminar  counterflow  model  with  the  shifting  equilibrium  model  indi¬ 
cates  that  differential  diffusion  effects  have  little  influence  on  the  Rayleigh/mixture  fraction  relationship 
These  results  indicate  that  mixture  fraction  measurements  can  lead  to  accurate  results,  both  for  highly  tur¬ 
bulent  diffusion  flames  and  laminar  strained  flamelets,  of  significant  value  to  the  modeling  of  these  flames 
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Captions 

Fig  l 

Fig  2 
Fig  3 
Fig  4 

Fig  5 
Fig  6 

Fig  7 

Fig8 

Fig  9 
Fig  10 

Fig  11 


Normalized  Rayleigh  signal  prediction  versus  mixture  fraction  for  a  fuel  of  SS^chelium 
,l5%hydrogen  burning  in  the  air,  calculated  using  a  shifting  equilibrium  model. 

Rayleigh  scattering  experimental  set  up  for  a  laminar  coflow  diffusion  flame. 

Normalized  Rayleigh  intensity  profile  obtained  at  10mm  from  the  burner  exit. 

Variance  of  the  Rayleigh  signal  versus  the  average  Rayleigh  intensity  signal  obtained  at  10mm 
of  the  laminar  burner  exit. 

Reaction  zone  structure  of  the  Tsuji  flame  for  a  strain  rate  a  =  10«-1. 

Normalized  Rayleigh  intensity  profiles  deduced  from  Eq.(l)  using  the  reaction  zone  prediction 
of  Fig.  5. 

Normalized  Rayleigh  intensity  profiles  comparisons  rescaled  using  the  maximum  gradient 
thickness. 

Comparison  between  predicted  and  measured  normalized  Rayleigh  intensity  profiles  in  mixture 
fraction  space. 

Reaction  zone  of  the  Tsuji  flame  for  a  strain  rate  a  ■*  190«-> 

Influence  of  strain  rate  on  the  normalized  Rayleigh  intensity  profile.  (1)  a  =*  10s-1  (2) 
«  -  190s-1. 

Comparison  of  the  normalized  Rayleigh  intensity  profiles  in  mixture  fraction  space  for  both 
laminar  counterflow  and  shifting  equilibrium  models. 
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ABSTRACT 

Detailed  scalar  and  velocity  measurements  have  been  made  in 
the  near  wake  region  of  bluff  body  stabilized,  lean,  premixed 
ethylene/air  flames  at  a  cold  flow  Reynolds  number  of  3800.  Two 
cases  were  studied  at  equivalence  ratios,  <b  =  0.63  and  <t>  —  0.54 
(close  to  the  lean  blow-off  limit).  Transverse  velocity  spectra  and 
velocity  fluctuation  patterns  characteristic  of  vortex  shedding  were 
observed  behind  the  rectangular  flame  holder  (3  mm  by  6  mm ). 
The  recirculation  zone  for  the  leaner,  less  stable,  flame  is  longer 
(65  mm  )  than  the  richer  case  (15  mm )  contrary  to  common  scaling 
assumptions.  This  is  probably  due  to  the  energetic,  low  frequency 
(<  70  Hz),  irregular  shedding  found  in  the  near  wake  of  the 
leaner  flame  which  may  be  an  indicator  of  impending  blow-off. 


1.  Introduction 

The  anchoring  of  premixed  flames  by  bluff  bodies  has  found 
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important  applications  in  the  stabilizing  of  high  speed  combustion 
in  ramjets  and  turbojet  afterburners.  It  has  also  been  suggested 
as  a  possible  mode  of  flame  holding  in  supersonic  combustion  dev¬ 
ices.  The  flame  holder  region  in  a  supersonic  combustor  is  likely 
to  be  subsonic  but  whether  currently  available  models  and  correla¬ 
tions  are  applicable  in  this  regime  is  not  known.  Furthermore,  in 
order  to  develop  flexible  numerical  modeling  strategies  more 
detailed  information  concerning  the  structure  of  the  stabilization 
region  is  essential. 

Much  of  the  early  experimental  work  on  flame  stabilization 
was  parametric  in  nature  and  directed  towards  the  determination 
of  blow-off  limits.  A  variety  of  intuitive  physical  models  were 
introduced  to  correlate  the  experimental  data.  Zukoski  and  Mar¬ 
ble  (1956),  for  example,  made  use  of  a  Damkohler  number  in  an 
ignition  delay  model  where  the  critical  approach  flow  velocity  at 
blow-off  is  reached  when  the  ignition  time  of  the  reactant  mixture 
exceeds  the  time  taken  for  a  fluid  element  to  flow  past  the  recircu¬ 
lation  zone.  Others  have  treated  the  recirculation  zone  as  a  well 
stirred  reactor  (Nicholson  and  Field,  1949)  or  compared  non- 
dimensional  heat  transfer  and  chemical  reaction  rates  (Spalding, 
1957).  These  studies  were  broadly  successful  in  generating  reason¬ 
ably  robust  correlations  but  provided  only  limited  insight  into  the 
detailed  structure  of  the  recirculation  zone  and  the  mechanism  of 
flame  stabilization.  The  most  significant  constraint  on  this  work 
was  the  range  of  diagnostic  techniques  available  which  were  res¬ 
tricted  largely  to  Schiieren  and  simple  scalar  measurements.  More 
recent  investigations  (Fuji  and  Eguchi,  1981)  using  point  measure¬ 
ment  techniques  present  a  more  complete  picture  but  the  charac¬ 
ter  of  the  flow  in  the  recirculation  and  flame  stabilization  region  is 
still  unclear.  There  is  little  evidence,  for  example,  of  vortex  shed¬ 
ding  in  bluff  body  stabilized  flames  although  the  Reynolds  number 
generally  found  in  hot  wakes  corresponds  to  the  vortex  street 
regime.  Cantwell  and  Coles  (1983),  for  instance,  have  observed 
vortex  streets  in  isothermal  flows  at  Reynolds  numbers  upto  10®. 

In  this  paper  the  results  of  a  detailed  experimental  study  of 
unconfined  premixed  turbulent  flame  stabilization  by  bluff  bodies 
are  reviewed  and  discussed.  Three  cases  are  investigated:  (a)  a 
‘well’  stabilized,  lean  flame,  (b)  a  flame  close  to  the  lean  blow-off 
limit  and,  for  purposes  of  comparison,  (c)  an  isothermal  flow.  The 
results  include  two-component  velocity  data  measured  by  laser 
Doppler  anemometry  (LDA),  gas  density  by  Rayleigh  scattering 
and  reaction  intensity  by  an  ionization  probe. 
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2.  Experimental  detail* 

A  schematic  of  the  burner  configuration  is  presented  in  figure  1. 
A  vertical  jet,  51  mm  in  diameter,  of  premixed  ethylene/air  is  sur¬ 
rounded  by  a  co-flowing,  co-axial  air  jet  which  shields  the  reactant 
stream  from  interaction  with  the  ambient  air.  The  flow  is 
unconfined  and  no  acoustic  resonances  are  observed  in  the 
approach  flow  spectrum.  5%  isotropic  turbulence  of  integral 
length  scale  2  mm  .  is  generated  by  a  square  grid  of  1  mm  cylindri¬ 
cal  rods  with  5  mm  spacing  placed  50  mm  upstream  of  the  nozzle 
exit.  The  flame  holder,  a  bar  of  6  4  mm  by  3.2  mm  cross  section, 
is  placed  75  mm  downstream  of  the  nozzle  exit  transverse  to  the 
mean  flow  direction.  Coordinate  directions,  measured  from  the 
top  center  of  the  flame  holder,  are  Z  in  the  axial  direction  ( U 
velocity)  and  X  in  the  transverse  direction  ( V  velocity). 

Flame  holding  at  two  equivalence  ratios  was  investigated,  a 
very  lean  flame,  close  to  lean  blow-off  at  equivalence  ratio  4>  — 
0.54  and  a  richer  flame,  <t>  =  0.63,  the  ‘well’  stabilized  case.  The 
Damkohler  numbers  for  both  of  the  flames  are  much  greater  than 
unity  indicating  that  the  flame  brush  away  from  the  stabilization 
region  may  be  characterized  as  a  wrinkled  laminar  flame  zone. 
Isothermal  measurements  were  also  made  and  experimental  condi¬ 
tions  are  referred  to  by  subscripts  /,  r  and  »  for  the  lean,  richer 
and  isothermal  cases  respectively.  The  freestream  velocity  is  6 
m  ft  in  all  cases  which  gives  a  cold  flow  Reynolds  number  of  3800 
based  on  the  wake  width  (defined  as  the  distance  between  n^, 
peaks). 

A  two-color  LDA  system  was  used  to  measure  the  axial  and 
transverse  velocity  components  and,  as  both  the  components 
undergo  flow  reversal  in  the  recirculation  zone  behind  the  flame 
stabilizer,  all  the  beams  were  frequency  shifted  to  give  a  12m  /» 
offset.  Velocity  bias  effects  were  corrected  by  weighting  the  data 
by  the  time  between  signal  validations.  Mean  and  rms  velocities 
were  calculated  from  8192  data  and  spectra  of  the  transverse  velo¬ 
city  component  were  obtained  bv  sampling  the  analog  signal  at  5 
kHz. 

Point  determinations  of  the  gas  density  were  obtained  by  Ray¬ 
leigh  scattering  with  a  probe  volume  of  0.1  mm  x  0.1  mm  .  The 
variation  in  scattering  cross-section  across  the  flame  front  is 
approximately  5%  in  lean  ethylene/air  flames  and  so  the  Rayleigh 
intensity  data  can  be  directly  interpreted  as  gas  density.  Mean 
density  statistics  and  scalar  spectra  were  calculated. 
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The  total  ion  concentration  across  the  flame  brush  which  is  a 
measure  of  the  reaction  intensity  is  determined  by  an  uncooled 
double  Langmuir-type  probe  with  an  active  length  of  5  mm  . 


3.  Results  and  Discussion 
3.1.  Mean  velocity  field 

Figure  2  a,b  show  the  right  hand  side  of  the  mean  velocity 
fields  for  the  two  flame  cases,  In  these  figures  the  magnitude  of 
the  vectors  may  be  best  appreciated  if  it  is  recalled  that  the  velo¬ 
city  in  the  freestream  is  6  m  jt  .  Also  marked  on  the  figures  are 
the  mean  positions  of  the  turbulent  flame  zone  which  are  derived 
from  Rayleigh  scattering  data.  The  most  immediately  observable 
difference  between  the  two  cases  is  the  length  of  the  recirculation 
zone,  L  .  For  the  leaner  case  L(  =  65  mm  which  is  twice  as  long 
as  the  richer  case,  Lr  =  32  mm  ,  indicating  that,  near  lean  blow- 
off,  this  length  is  not  a  good  indicator  of  flame  stability.  Lt  and 
Lf  are  both  longer  than  the  isothermal  case  which  is  in  agreement 
with  other  studies  (Chen  et  al,  1988). 

The  point  at  which  the  flame  can  be  considered  fully  ignited 
can  be  deduced  from  the  position  of  the  flame  zone  relative  to  the 
shear  layer  between  the  freestream  and  the  recirculation  zone.  In 
the  richer  case  the  flame  is  able  to  propagate  out  of  the  shear 
layer  at  0.5Lr  but  in  the  lean  case  this  takes  place  downstream  of 
the  rear  stagnation  point  and  so  the  convection  times  to  self  pro¬ 
pagation  differ  by  a  factor  of  five.  The  chemical  times,  estimated 
from  the  laminar  flame  thickness  and  laminar  burning  velocity, 
however,  differ  by  a  factor  of  three  indicating  that  problems  may 
be  encountered  if  stability  criteria  are  based  on  comparisons  of 
these  times.  The  establishment  of  a  propagating  flame  is  also 
observable  in  the  affect  of  the  flame  on  the  mean  velocity  vectors. 
Downstream  of  the  recirculation  zone  and  away  from  the  wake 
region  the  velocity  vectors  are  deflected  towards  the  centerline  as 
the  reactants  are  burned,  expand  and  so  are  accelerated  at  the 
flame  front.  Since  the  included  angles  for  the  flames  are  small  (0, 
=  7 *,  9,  =  16*)  this  acceleration  is  primarily  in  the  transverse 
direction.  Where  the  flame  is  in  the  shear  layer,  however,  and  not 
fully  established  little  of  no  deflection  occurs:  this  is  clearly  seen  in 
the  leaner  case.  It  may  also  be  noted  that  the  growth  of  in  the 
flame  brush  thickness,  6? ,  is  small  in  the  leaner  case  indicating 
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little  interaction  with  the  freestream  turbulence  which,  in  a  freely 
propagating  flame,  causes  the  flame  brush  to  grow  by  turbulent 
diffusion  (Goix  et  al,  1989)  as  can  be  seen  in  the  richer  case. 
There  is  an  increase  of  25%  in  (6T)t  between  Z  —  20  mm  and  Z 
=  40  mm  while  over  the  same  distance  (<5r)f  increases  by  a  factor 
of  3. 


3.2.  Velocity  Fluctuations  and  Spectra 

Velocity  measurements  made  at  the  nozzle  exit  show  that  the 
turbulence  structure  is  isotropic  and  has  an  energy  spectrum 
characteristic  of  grid  turbulence.  The  structure  of  the  recircula¬ 
tion  zone  in  isothermal  conditions  will  be  analysed  briefly  as  it  will 
aid  in  the  interpretation  of  the  velocity  fields  of  the  combusting 
cases.  The  cold  flow  Reynolds  number  of  3800  indicates  that  vor¬ 
tex  shedding  behind  the  bar  is  to  be  expected  and  the  transverse 
velocity  spectrum  has  a  prominent,  narrow  peak  at  100  Hz.  The 
Strouhal  number  of  0.16  is  lower  than  that  observed  for  cylindrical 
bluff  bodies  (0.2)  but  this  may  be  expected  due  to  the  rectangular 
shape  of  the  bar.  The  energy  of  the  peak  in  the  spectra  has  a  max¬ 
imum  at  the  end  of  the  recirculation  zone,  Z  =  25  mm  ,  and 
thereafter  decays  slowly.  The  vortex  shedding  also  has  a  marked 
effect  on  the  distribution  of  the  velocity  fluctuations  and  figures  3 
a,b  show  contour  maps  of  the  axial  and  transverse  velocity 
fluctuations  for  the  isothermal  case.  The  maximum  axial  fluctua¬ 
tions,  figure  3  a,  are  found  in  the  shear  layer  between  the  recircu¬ 
lation  zone  and  the  freestream.  The  maximum  transverse  fluctua¬ 
tions,  figure  3  b,  however,  occur  at  the  centerline  and  are  due  to 
the  alternately  shed  vortices  being  swept  across  the  centerline. 
These  phenomena  are  characteristic  of  vortex  shedding  behind 
bluff  bodies  (Cantwell  and  Coles,  1983)  and  the  results  are  sum¬ 
marized  in  Table  I  where  L  is  the  recirculation  zone  length  and  / 
the  shedding  frequency. 
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Table  I:  Recirculation  Zone  Data 

4> 

Re 

L 

/ 

«mu 

mm 

Hz 

m  ft 

m  ft 

0.0 

3800 

27 

100 

1.94 

2.23 

0.54 

150 

65 

70 

1.78 

1.00 

0.63 

110 

32 

80 

1.00 

0.64 

Figures  4  a,b  and  5  a,b  compare  the  axial  and  transverse  velo¬ 
city  fluctuation  contour  plots  for  the  two  flame  cases.  The  center 
of  the  flame  brush  is  also  marked  on  these  figures.  In  considering 
these  results  it  is  important  to  distinguish  the  various  causes  of 
velocity  fluctuations  above  those  found  in  the  freestream  (approxi¬ 
mately  0.2  m  /» ).  In  regions  where  the  instantaneous  flame  front 
is  flapping  across  the  measurement  point  the  velocity  difference 
between  the  burnt  and  unburnt  gases  will  produce  an  apparent 
increase  in  the  velocity  fluctuations  (Cheng,  1984).  As  mentioned 
above,  this  acceleration  occurs  predominately  in  the  transverse 
direction  and  its  effect  can  be  seen  most  clearly  ;n  the  v^,  contour 
plot  for  the  richer  case,  figure  4  b,  where  significant  fluctuations 
along  the  mean  flame  position  start  to  appear  as  the  flame 
separates  from  the  shear  layer  at  Z  =15  mm  .  Axial  fluctuations 
in  the  richer  case,  figure  4  a,  only  increase  in  the  flame  zone  where 
the  m-an  flame  angle  has  changed  sufficiently  for  the  acceleration 
to  contribute  to  the  axial  velocity  component.  The  leaner  flame, 
figures  5  a,b,  is  markedly  different  and  the  presence  of  the  flame 
has  little  effect  on  the  velocity  fluctuations  in  either  component. 
This  underlines  the  largely  passive  nature  of  the  flame  during  the 
initial  phase  of  its  development. 

The  other  major  region  of  velocity  fluctuation  increase  in  both 
the  combusting  cases  is  the  recirculation  zone.  The  maximum 
velocity  fluctuations  in  the  recirculation  zone  are  given  in  Table  I. 
The  maxima  are  somewhat  reduced  from  the  isothermal  values 
and  this  is  probably  due  to  the  volumetric  expansion  and  the 
much  reduced  Reynolds  number.  Setting  aside  those  factors  asso¬ 
ciated  with  the  flame  zone  which  were  considered  above,  it  is  clear 
that  the  distributions  of  velocity  fluctuations  are  very  similar  to 


the  isothermal  case.  The  position  of  the  peak  velocity  fluctuations 
and  its  evolution  occur  with  the  same  pattern  and  indicate  the 
existence  of  vortex  shedding  in  the  near  wake  region.  Figure  6 
shows  the  transverse  velocity  spectra  on  the  centerline  for  the  two 
flames  at  Z  =  20  mm  .  The  spectrum  of  the  richer  flame  at  the 
centerline  at  Z  =20  mm  has  a  sharp,  energetic  peak  at  80  Hz 
which  gives  a  Strouhal  number  of  0.1.  An  isothermal  value  (Oka- 
jima,  1982)  for  square  cylinders  of  0.13  has  been  measured  but 
estimation  of  the  Reynolds  number  in  the  wake  is  not  straightfor¬ 
ward.  The  bluff  body  is  not  in  uniform  flow  but,  using  the 
kinematic  viscosity  in  the  hot  wake,  the  width  of  the  flame  holder 
and  the  freestream  velocity,  a  value  of  110  is  obtained.  These 
results,  with  the  distribution  of  the  velocity  fluctuations,  indicate 
the  existence  of  vortex  shedding  behind  the  flame  holder.  The 
spectrum  of  the  leaner  flame,  which  has  a  slightly  higher  Reynolds 
number  of  150.  has  a  small  peak  at  70  Hz.  The  main  feature  of 
the  spectrum,  however,  is  the  large  amount  of  energy  below  the 
shedding  frequency  which  shows  that  the  shedding  has  become 
highly  irregular.  This  may  explain  the  large  size  of  the  recircula¬ 
tion  zone  and  may  be  viewed  as  a  precursor  to  blow-off.  It  is  of 
interest  to  compare  the  spectrum  of  the  density  fluctuations  at  the 
same  height  above  the  flame  holder  but  at  the  center  of  the  flame 
brush,  figure  7.  The  large  velocity  fluctuations  do  not  strongly 
affect  the  flame  oscillations  and  so  a  measurement  technique  such 
as  Schlieren,  which  detects  density  gradients,  would  not  reveal  the 
unsteady  nature  of  the  recirculation  zone. 

3.3.  Scalar  measurements 

The  Rayleigh  measurements  have  been  used  to  delineate  the 
xtent  of  the  flame  zone  and  to  derive  spectra  of  density  fluctua¬ 
tions.  Further  insight  into  the  development  of  the  flames  may  be 
obtained  from  a  consideration  of  the  probability  density  functions 
(pdf)  of  the  Rayleigh  scattering  intensity.  Figure  8  a,b  shows  the 
density  pdfs  at  the  center  of  the  flame  brush  for  the  two  flames 
from  Z  —  5  mm  to  Z  —  50  mm  .  When  the  flame  brush  is  large 
compared  to  the  laminar  flame  thickness,  ie  in  the  wrinkled  lam¬ 
inar  flame  regime,  the  pdf  will  be  strongly  bimodal  with  a  low  pro¬ 
bability  of  detecting  densities  other  than  those  characteristic  of 
the  burned  or  unburned  gases.  This  can  be  clearly  seen  for  posi¬ 
tions  Z  >  15  mm  of  the  richer  case  when  the  flame  has  pro¬ 
pagated  out  of  the  shear  layer.  The  leaner  flame  remains  largely 
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within  the  shear  layer,  the  bimodality  is  much  less  apparent,  and 
close  to  the  flame  holder  a  broad  distribution  of  densities  is 
observable.  At  this  position  ( 6r)(  =1.5  mm  which  is  of  the  same 
order  as  the  laminar  flame  thickness  (0.8  mm  )  and  therefore  this 
distribution  could  arise  from  a  weakly  fluctuating  laminar  flame 
The  ion  probe  results,  figure  9,  show  that  this  is  not  the  case  and 
that  very  little  chemical  reaction  is  taking  place  at  these  positions. 
It  is  of  interest  to  note  that  the  integrated  ion  current,  in  the 
richer  case,  rises  to  a  maximum  at  the  point  where  the  flame 
leaves  the  shear  layer. 


4.  Conclusions 

The  understanding  of  bluff  body  flame  holding  has  largely 
been  obtained  by  the  use  of  dimensional  reasoning  based  on  rela¬ 
tively  simple  physical  models  which  assume  steady  recirculation 
behind  the  flame  stabilizer.  Whether  these  concepts  are  applicable 
to  flame  anchoring  in  super  sonic  combustors  is  unclear. 

A  detailed  experimental  investigation  of  lean  ethylene/air 
flame  stabilization  has  been  presented  here  as  a  step  towards  a 
more  complete  understanding  of  flame  holding.  This  study  shows 
that  vortex  shedding  which  was  observed  in  the  near  wake  o.  the 
flame  holder  can  have  significant  affects  on  the  flame  stability 
Close  to  the  lean  blow-off  limit  the  shedding  becomes  highly  irreg¬ 
ular  and  is  probably  the  cause  of  the  large  increase  in  the  recircu¬ 
lation  zone  length  and  the  higher  velocity  fluctuations.  These 
phenomena  may  be  indicators  of  impending  blow-off  and  may 
indeed  be  responsible  for  it. 
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Figure  1.  Burner  Configuration. 


Figure  2.  Mean  Velocity  Field. 
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Figure  3.  Velocity  Fluctuations  a)  axial  b)  transverse 
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.  Figure  4.  Velocity  Fluctuations  a)  axial 
b)  transverse. 
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Figure  5.  Velocity  Flucuations  a)  Axial 
b)  transverse. 
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Figure  8. 
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Experimental  Study  of  Turbulent  Counterflow  Diffusion 
Flame  Structure 
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ABSTRACT 

A  Hydrogen-Helium  mixture  has  been  chosen  to  investigate  the  dilution  effects  on  the  structure  of  the  diffu¬ 
sion  flame.  To  investigate  the  combustion  effects  on  the  turbulent  mixing,  reacting  and  non  reacting  conditions  have 
been  studied  at  the  same  Reynolds  number.  In  order  to  study  the  reaction  zone  structure,  high  speed  tomography 
based  on  Mie  scattering  has  been  employed  using  a  copper  vapor  laser  and  a  Fastax  high  speed  camera.  Different 
seeding  techniques  have  been  used  to  visualize  both  the  turbulent  air  and  fuel  jets.  The  tomographic  records  were 
digitized  and  recorded  in  a  digital  computer  for  statistical  treatment  Significant  differences  in  the  wrinkle  scales 
between  the  reacting  and  the  non  reacting  flows  were  found.  A  fractal  statistical  analysis  of  the  tomography  records 
has  been  done  to  quantify  these  differences.  Seeding  of  both  fuel  and  air  jets  provided  a  mean  for  the  evaluation  of 
the  reaction  zone  thickness.  From  the  time  resolved  tomographic  records  an  evaluation  of  the  stretch  of  the  reaction 
zone  boundaries  has  been  performed.  Local  flame  extinction  and  reignition  have  been  observed  for  different 
H2/Helium  fuel  mixture.  Further  statistic  on  reaction  zone  thickness  are  still  in  progress. 


291 


(4 

% 


P 1 8  3  .  Direct  Measurement  of  Mixture  F raction  in 

Reacting  Flow  Using  Rayleigh  Scattering 
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Abstract 

The  mixture  fraction  variable,  |,  is  very  useful  in  describing  reaction  zone  structure  in  non-premixed  flames. 
Extinction  limits  and  turbulent  mixing  are  often  described  as  a  function  of  this  variable.  Experimental  evaluation  of 
^  is  valuable  in  understanding  the  influence  of  turbulent  mixing  on  the  chemistry  process.  Historically,  the  evalua¬ 
tion  of  mixture  fraction  in  combusting  flow  has  required  multiple  concentration  measurements.  Here,  a  fuel  has 
been  designed  to  permit  the  measurement  of  the  mixture  fraction  using  only  Rayleigh  scattering.  A  Rayleigh  inten¬ 
sity/  mixture  fraction  correspondence  has  been  obtained  experimentally  in  a  laminar  coflow  flame.  The  influence  of 
strain  rate  and  differential  diffusion  effects  has  been  investigated  using  a  laminar  counterflow  diffusion  flame  and 
shifting  equilibrium  models.  The  results  obtained  from  these  comparisons  are  very  encouraging  and  suggest  that  the 
Rayleigh/  mixture  fraction  correspondence  established,  is  valid  under  both  turbulent  mixing  and  laminar  strained 
flamelet  combustion  regimes. 
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